Ab-initio computation of neutron-rich oxygen isotopes 
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We compute the binding energy of neutron-rich oxygen isotopes and employ the coupled-cluster 
method and chiral nucleon-nucleon interactions at next-to-next-to-next-to-leading order with two 
different cutoffs. We obtain rather well-converged results in model spaces consisting of up to 21 
oscillator shells. For interactions with a momentum cutoff of 500 MeV, we find that 28 is stable 
with respect to 24 0, while calculations with a momentum cutoff of 600 MeV result in a slightly 
unbound 28 0. The theoretical error estimates due to the omission of the three-nucleon forces and 
the truncation of excitations beyond three-particle-three-hole clusters indicate that the stability of 
28 cannot be ruled out from ab-initio calculations, and that three-nucleon forces and continuum 
effects play the dominant role in deciding this question. 
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Introduction. The neutron drip line marks the limits of 
stability of neutron-rich isotopes. At present, this line is 
well established only in the lightest elements, as the cross 
section for the production of extremely neutron-rich nu- 
clei decreases dramatically as one moves away from nu- 
clei in the valley of beta stability (for a recent review 
see for example Ref. jlj)- At present, 24 is the "last" 
known stable neutron-rich oxygen (Z = 8) isotope, and 
25 O is known to decay under the emission of one neutron 
0. The "next" neutron-rich oxygen isotope 26 O has not 
been observed experimentally [3|, [j] , and systematics for 
its production cross section suggest that it should have 
been seen if it were a stable nucleus. Similar estimates 
suggest that the isotope 28 is unstable [5J]. Thus, exper- 
iment puts the neutron drip line at 24 O. This is remark- 
able since 31 F is the most neutron-rich fluorine (Z = 9) 
isotope [f|. Thus, the addition of a single proton appar- 
ently shifts the drip line by six neutrons. 

The theoretical determination of the neutron drip line 
is a challenging task as well. Near the neutron drip 
line, small uncertainties in the nuclear interaction are 
enhanced due to the extreme isospin and the proximity 
of the continuum. Several theoretical studies have ad- 
dressed the structure of neutron-rich oxygen and fluorine 
isotopes. The employed methods and theoretical pre- 
dictions differ considerably. The sd-shell model, based 
on the Brown- Wildenthal USD interaction and the fi- 
nite range droplet model, predicts that 26 is stable Q- 
Within the sd-pf shell model, the present experimen- 
tal situation of an unstable 26 O can be reproduced af- 
ter a modification of the interaction Q. Within this 
model, 28 is unbound by about 1 MeV. Within the 
same model space, but a different interaction, Poves and 
Retamosa obtained a stable 31 F and a stable 28 0. 
Shell- model descriptions of neutron-rich oxygen isotopes, 
including the coupling to the continuum, were given in 
Refs. 13, U, 13] • Within the latter approach 12 1, two 



slightly different phenomenological sc?-shell interactions 
are employed for oxygen isotopes close to and far away 
from the valley of beta stability, respectively. This leads 
to the result that 26 is unstable with respect to 24 0, 
while 28 O is unstable with respect to the emission of two 
and four neutrons. Clearly, the present theoretical sit- 
uation does not have the desired predictive power, and 
calculations suffer from uncertainties in the knowledge 
of the interaction and from the difficulty to quantify how 
these uncertainties propagate in the quantum many-body 
problem. This is an opportunity for ab-initio calculations 
to address these challenges. 

Ab-initio calculations have been very successful in light 
nuclei [13|, uA, ll5j, Il6l Il7l l , a nd have recently also been 
extended to unbound |l8l . ItqI |20| and medium-mass iso- 
topes plj . In this paper we present ab-initio calcula- 
tions for the neutron-rich oxygen isotopes 22 > 24 > 28 0, and 
employ nucleon-nucleon interactions from chiral effective 
field theory (EFT) |1 [H, 0, [H, HI- These interac- 



tions are rooted in quantum-chromo-dynamics and in- 
clude pion exchange and zero-range forces. The power 
counting, i.e., the systematic expansion of the interac- 
tion in terms of ratios of the probed momentum scale 
Q over the cutoff A v is an important asset. In finite 
nuclei, Q is about 



x is an important asset. 
200 MeV [2l|, while the cutoffs we 



employ are A x = 500 MeV and A x = 600 MeV, respec- 
tively. The variation of our results with the cutoff allows 
us to quantify uncertainties that are due to the omis- 
sion of (short-ranged) three-nucleon forces. We employ 



the coupled-cluster method |28j, |29|, |3fj, \M, lM |33j for 
the solution of the quantum many-body problem. This 
method scales gently with the system size, and can accu- 
rately compute the binding energies of nuclei with closed 
subshells. In particular, the possibility to employ large 
model spaces avoids the need for a secondary renormal- 
ization of the chiral interactions from EFT for nuclei such 



as oxygen and calcium isotopes [21] 
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This paper is organized as follows. We briefly intro- 
duce the interactions and methods we employ and then 
present the results of our calculations. 

Interaction, model space and coupled- cluster method: 
We employ the chiral nucleon-nucleon interaction by 
Entem and Machleidt [25| at next-to-next-to-next-to- 
leading order (N 3 LO). This includes terms up to order 
(Q/A x ) 4 in the power counting of the nucleon-nucleon 
interaction. The interaction has a high-momentum cut- 
off of A x = 500 MeV, and a version with cutoff A x = 
600 MeV is also available [26]. The low-energy con- 
stants of the chiral potentials were determined by fits to 
the two-nucleon system. We neglect three-nucleon forces 
that already appear at next-to-next-to-leading order and 
thereby introduce uncertainties of the order (Q/A x ) 3 . 
As physics must be independent of the cutoff (or renor- 
malization scale), any cutoff-dependence in our results 
quantifies the uncertainty due to omitted contributions 
of short-ranged three- nucleoli forces and forces of higher 
rank. The intrinsic Hamiltonian is 



H = T -Tr: 



V(A X ) 



Here T, f cm , and V(A X ) denote the kinetic energy of 
the A-body system, the kinetic energy of the center-of- 
mass coordinate, and the chiral nucleon-nucleon inter- 
action with momentum cutoff A x , respectively. The in- 
trinsic Hamiltonian ([T]) is translationally invariant and 
does not depend on the center-of-mass coordinate. We 
express the Hamiltonian in a single-particle basis of the 
spherical harmonic oscillator. Our model-space param- 
eters are the oscillator spacing hw of our single-particle 
basis, and the maximal excitation energy (TV + 3/2)Huj 
of a single-particle state, that is the number of complete 
oscillator shells is N + 1. As a first step towards the so- 
lution of the many-body problem, we solve the spherical 
Hartree-Fock equations and transform the Hamiltonian 
to this basis. In the second step, we employ the coupled- 
cluster method. In drip-line nuclei, the outermost nu- 
cleons move in orbitals close to the scattering threshold, 
making the nuclear wave function exhibit halo-like struc- 
tures and sometimes even ground states embedded in the 
continuum. The presence of the scattering continuum in 
such exotic nuclei makes the use of the oscillator basis 
not ideal. The unrealistic Gaussian falloff of the oscil- 
lator wave functions makes convergence slow for nuclei 
with dilute matter distributions. However, the Gamow- 
Hartree-Fock (see, e.g., Ref. [Hj]) yields occupied single- 
particle states with nonphysical positive imaginary parts. 
This difficulty is due to the relatively "hard" interaction 
we employ. To avoid this problem, we choose to stay 
within the oscillator basis, but employ very large model 
spaces for an improved description of the tails of the ra- 
dial wave function. 

The nuclear many-body problem is solved with the 
coupled-cluster method [II, ^ 0, HH, [H [H, [H, [H, H3, 



1 3 81 ] . This approach is based on the similarity transforma- 
tion of the normal-ordered intrinsic Hamiltonian Hn , 
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Here, the Hamiltonian is normal-ordered with respect to 
a product state which serves as a reference. The 
particle-hole cluster operator 



T = T 1 +T 2 +T 3 + 



(2) 



is defined with respect to this reference state. It is a sum 
of the fc-particle fc-hole (fep-feh) cluster operators 



T, 



1 



(fc!) 



H2 



E 



xoi...o fe -t 
b i!...t k u ai 



Q'a.k^ik ■ ■ - a i 



(3) 

We use the convention that i,j,k,... label occupied 
single-particle orbitals while a, b, c, . . . label unoccupied 
orbitals. We truncate the cluster operator beyond the T 2 
level and employ ACCSD(T) [39|,[40] as an approximation 
for the T3 clusters. The unknown cluster amplitudes tf 
and t£ in Eq. © are determined from the solution of 
the coupled-cluster equations 



(4) 



Here \ipf) = a+aj| 



and \ipfj) = a a a t a jO,i\ip) ar e lp-lh 



and 2p-2h excitation of the reference state, respectively. 

The nonlinear coupled-cluster equations ^ are solved 
iteratively, and the correlation energy of the ground state 
is computed from 



A^ccsd = bl)\H\ 



(•5) 



We employ a spherical formulation of coupled-cluster 
theory where the cluster operator T is a scalar under 
This formulation reduces considerably the 
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number of unknowns and permits us to explore model 
spaces exceeding 20 major oscillator shells. 

Let us briefly summarize the essential properties of the 
coupled-cluster method. First, the method fulfills Gold- 
stone's linked cluster theorem and therefor yields size- 
extensive results, i.e., the error due to the truncation is 
linear in the mass number A. Size extensivity is an im- 
portant issue when approximate solutions to all but the 
lightest nuclei are sought (3^, EJ. Second, the compu- 
tational effort scales gently (i.e., polynomial) with the 
system size. The method has met benchmarks in light 
nuclei [12, E^]. We neglect three-nucleon forces since 
their application within the coupled-cluster method is 
presently limited to smaller model spaces [44(. 

For a more precise computation of the correlation en- 
ergy, we consider corrections due to triples excitations T3 
within the ACCSD(T) approximation. For this purpose, 
we solve the left eigenvalue problem 



(V>|AfT = E{ij}\K 



(6) 
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of the similarity-transformed Hamiltonian H. Here, A 
denotes the de-excitation cluster operator 



A = 1 + Ai + A 2 



(7) 



and 



Ai = A « a ° a l ' ^ = jJ2 ^%a b a a a\a) ■ (8) 



i,j,a,b 



The unknowns and X^ b result from the ground-state 



solution of the left eigenvalue problem (|6J). They are 
utilized together with the cluster amplitudes tf and tfj 
to compute the energy correction due to triples clusters 
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Here, Fh p denotes the part of the normal-ordered one- 
body Hamiltonian that annihilates particles and creates 
holes, while e^r = f u + f j} + f kk - f aa - f bb - f cc is 
expressed in terms of the diagonal matrix elements of the 
normal-ordered one-body Hamiltonian F. The subscript 
C denotes the connected part of the operator, and IV'yfc) 
is a 3p-3h excitation of the reference state. 

Results: We consider the nuclei 16 > 22 > 24 > 28 and com- 
pute their ground-state energies within the ACCSD(T) 
approximation for chiral interactions with cutoffs of A x = 
500 MeV and A x = 600 MeV, respectively. FiguresCDand 
[2] show the results as a function of the oscillator spacing 
Huj of the single-particle basis, and parametrized by the 
number of major oscillator shells N + 1 for 24 O, and 28 O 
with two chiral cutoffs A x , respectively. Note that the 
results are reasonably well converged with respect to the 
size of the model space. Note also that the "harder" inter- 
action with cutoff A x = 600 MeV requires a larger model 
space to reach an acceptable convergence. The results for 
16 O and 22 O are of similar quality. Let us also comment 
on the separation of the center-of-mass coordinate and 
the intrinsic coordinates. Very recently, Hagen, Papen- 
brock, and Dean demonstrated that the coupled-cluster 
wave function is approximately a product of a transla- 
tionally invariant wave function and a Gaussian for the 
center-of-mass coordinate [45j. Thus, we do not worry 
about spurious contributions to the coupled-cluster wave 
function. 

Let us estimate the precision of our results. There are 
three sources of systematic errors, namely the truncation 
level of the coupled-cluster method, the finite size of the 
model space, and the error due to omitted contributions 
in the interaction. First, within the ACCSD(T) approxi- 
mation, three-particle-three-hole clusters are treated ap- 
proximately, and all excitation clusters of higher rank 
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FIG. 1: (Color online) Binding energy (within ACCSD(T)) 
for 24 from a chiral NN potential at order N 3 LO with high- 
momentum cutoffs A x = 500 MeV as a function of the oscil- 
lator spacing hui and the size of the model space. 
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FIG. 2: (Color online) Same as Fig. [T] except for 28 O and a 



chiral cutoff A x = 600 MeV. 



are neglected. Table Q] shows the different contribu- 
tions to the binding energy of neutron-rich oxygen iso- 
topes. Comparison of the CCSD correlation energy and 
the energy due to triples corrections shows that the lat- 
ter account for 10% (13%) of the former at a cutoff 
A x = 500 MeV (A x = 600 MeV). These ratios are 
found in similar coupled-cluster calculations of atoms and 
molecules, and experience in quantum chemistry (see for 
example Ref. (33|) suggests that the truncation of the 
cluster amplitudes beyond the triples corrections intro- 
duces an error of a few percent corresponding to an un- 
certainty of approximately 5 MeV. Second, we cannot 
treat an infinite model space, and (as shown in Figs. [1] 
and[2|), the convergence with respect to an increased size 
of the model space is at the level of a couple of MeV. 
Thus, the convergence with respect to the size of the 
model space introduces an error that does not exceed 
error estimates due to the truncation of the cluster am- 
plitudes. Third, by far the largest uncertainty is due to 
omissions in the nuclear interaction, as can be seen from 
a comparison of the results obtained with two different 
cutoffs. This uncertainty is of the order of 10-20 MeV, 
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and increases with increasing mass number. Note that 
the deviation from the experimental results is consistent 
with our error estimates. Overall, we are missing bind- 
ing energy compared to experiment. Thus, the net effect 
of the three-nucleon force is expected to be attractive 
for both cutoffs. Table |TT] gives the root-mean-square 
point matter radii using the chiral interaction with high- 
momentum cutoff A x = 500MeV. 

Radii were calculated using the Helmann-Feynman 
theorem within the ACCSD(T) approximation in 19 ma- 
jor oscillator shells. We compare with the effective point 
matter radii extracted from interaction cross sections us- 
ing the Glauber model in the optical limit approximation 
34j . Our calculated matter radii are smaller than those 
extracted from experiment. In our calculated radii we 
estimate an uncertainty at the order of ~ 0.1 fm from 
the model-space dependence. 

Let us also check whether our error estimates are 
consistent with the power-counting estimates from chi- 
ral EFT. Nogga confirmed that these estimates hold in 
light nuclei [461 1 . The omitted three-nucleon forces are 
of the order (V)(Q/A X ) 3 where Q is the typical mo- 
mentum scale and (V) is the expectation value of the 
two-body interaction. For nuclei in this mass region and 
a cutoff A x = 500 MeV, we have Q « 200 MeV and 
(V) w 33 ± 3 MeV per nucleon (taken from the ex- 
pectation values of the kinetic and potential energies in 
Ref. [2l|, respectively). This puts power-counting esti- 
mates from chiral EFT at about 2 MeV per nucleon, 
and our results are well within this estimate. While the 
absolute uncertainty on the binding energy is thus con- 
siderable, the differences in the binding energies of the 
considered isotopes (at fixed chiral cutoff A x ) is much 
closer to the experimental result. 

Let us focus on the binding energy of 28 O with respect 
to 24 O. While it would certainly be interesting to include 
26 O in this comparison, we cannot address this nucleus 
within the spherical coupled-cluster method due to its 
open-shell character. Recall, however, the experimental 
evidence 0, 3| against the stability of 26 O. This makes 
the comparison of 28 O and the last known stable isotope 
24 particularly interesting. Figure [3] shows that the 
ground-state energies relative to 22 O change little as one 
goes from 24 to 28 0. This is a remarkable result of our 
ab-initio calculations. In shell-model calculations with 
an 16 O core, the ground-state energies typically increase 
strongly (in absolute value) as more neutrons are added, 
and an adjustment of the interaction is necessary [l|. 
We study this phenomenon further by employing a sim- 
ilarity renormalization group (SRG) transformation (47| 
of the Hamiltonian with cutoff A x = 500 MeV. As we 
lower the smooth SRG momentum cutoff from 4.1 fm -1 
to 3.5 fm -1 , we find that the ground-state energy of 28 
decreases further relative to 24 0. Thus, a softening of 
the nucleon-nucleon interaction has to be compensated 
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46.33 
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-119.04 
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-168.49 
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AE :i 


-14.95 


-20.71 


-22.49 


-22.86 


E 


-111.91 


-130.89 


-138.04 


-136.71 


Experiment 


-127.62 


-162.03 


-168.38 





TABLE I: Contributions to the binding energy E (in MeV) 
in neutron-rich oxygen isotopes from chiral interactions with 
high- momentum cutoff A x . The contributions Eo, ASccsd, 
and AE3 denote the Hartree-Fock energy, the correlation en- 
ergy within the CCSD approximation, and the energy due 
to the employed triples correction, respectively. For le O and 
A x = 500MeV the results were taken in 19 major oscillator 
shells and at the energy minimum hu = 40MeV. For all other 
cases the results were obtained in the largest model spaces at 
fixed huj = 28 MeV. 
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2.54(2) 


2.88(6) 


3.19(13) 





TABLE II: Root-mean-square point matter radii (in fm) for 
neutron-rich oxygen isotopes from the chiral interaction with 
high-momentum cutoff A x = 500MeV. Oscillator frequencies 
as in Table [T] Experimental data from Ref. [33 ], 



by a three-nucleon force that yields less attraction (or 
even repulsion) in 28 O than in 24 O. There is no cutoff 
in this range that simultaneously would reproduce the 
experimental binding of 22 O and 24 O. 
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FIG. 3: (Color online) Ground-state energies of neutron-rich 

isotopes A rela 
two different cutoffs A x 



oxygen isotopes A relative to 22 O for chiral interactions with 



At a cutoff A v = 500 MeV we find that 28 O is bound 
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by about 2.7 MeV with respect to 24 0. However, the sit- 
uation is reversed at the higher cutoff A x = 600 MeV, 
and the difference is about -f.3 MeV. Given the uncer- 
tainties of our calculation as discussed in the preceding 
paragraph, it is presently not possible to reach a con- 
clusion regarding the existence of 28 O. However, the dis- 
cussion presented in the previous paragraph also makes 
clear that - within interactions from chiral EFT - the 
stability of 28 O depends mainly on the contributions of 
the three-nucleon force, and that even small contribu- 
tions can tip the balance in either direction. This is the 
main result of this paper. Our ab-initio calculations also 
suggest that the recent results from phenomenological 
shell-model approaches regarding the unbound character 
of 28 O might be viewed with caution. The combination 
of three-nucleon forces, the proximity of the continuum 
and the isospin dependence are a challenge for reliable 
theoretical predictions. 

In summary, we performed ab-initio calculations for 
neutron-rich oxygen isotopes employing chiral nucleon- 
nucleon interactions at order N 3 LO. We probed the ef- 
fects of missing physics (such as three-nucleon forces) by 
studying the cutoff dependence of our results, and es- 
timated the uncertainties due to the finite size of the 
model space and the truncation of the cluster operator. 
Our results show that the absolute binding energies have 
considerable uncertainties. However, the differences in 
binding energies are much closer to experiment. We find 
a small difference in the binding energies of 24 O and 28 O. 
Thus, our results cannot rule out a stable 28 with re- 
spect to 24 O. The cutoff-dependence of the results shows 
that three-nucleon forces are the dominant contributions 
that tip the balance. 
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